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Abstract 
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1 Introduction 



A number of differential equations providing mean-field descriptions for the ki- 
netics of particle aggregation have been the focus of much mathematical effort 
(see |21| for a review of recent mathematical developments.) An important spe- 
cial case of this type of equations are the well known coagulation equations first 
studied by the Polish physicist Marian Smoluchowski (1872-1917) I3U]: denot- 
ing by Cj = Cj(t) the concentration at time t of an agglomerate of j identical 
particles, and assuming binary aggregation following mass action law is the only 
process taking place, we obtain the discrete version of Smoluchowski's system 

i j — 1 oo 
Cj = ^ X] a 3-kMCj- k C k - Cj ^ a J,kCk (1) 
k=l k=l 

where the first sum in the right-hand side is defined to be zero when j = 1. The 
kinetic coefficients a^h measure the efficiency of the reaction between j— clusters 
and k— clusters to produce (j + fc)-clusters, and, as such, they should satisfy 
the minimal requirements of symmetry and nonnegativity a^fc = a^.j > 0. 
The precise form of the coefficients depends on the physical phenomena being 
modelled (see, eg Table 1].) 

A central problem in the study of the long-time behaviour of solutions to 
JQl is their convergence to some similarity profile ip, when j and t — > +oo, 

Cj (t) ~ S (t)-?iP(js(t)- a ) 

where a and (3 are appropriate positive exponents, and s(t) a scaling function. 

This is conjectured to occur, for a large class of initial conditions, in the 
case of homogeneous kinetic coefficients, i.e., those satisfying cij.k = K(j,k), 
with kernels K for which K(ux, uy) = u x K(x, y). For a very recent survey see 

E3 

From a rigorous point of view not much is known, although some results 
were previously available I18| . and a number of significant advances have 
very recently been made ^3 EJ EH1 ES] ■ 

In the present paper we are interested in studying this type of behaviour for 
a class of non-homogeneous coefficients that are also relevant in the applications: 
the Becker-Doring type coefficients satisfying Oj,*, = if min{j, k} > 1. With this 
type of restriction system Q is sometimes called the addition model [3 llfil , 
and is written in the form 



ci = — dicf — c±y ] a j c j 

3 = 1 

c 3 = Qj-iCiCj-l — Oj'CiCj, j>2 



(2) 



where a\ — and dj = 07,1 if j > 2. It is a special case of the Becker-Doring 

coagulation equations 0] (without fragmentation of clusters.) Physically this 
corresponds to cases where the only effective reactions are those involving at 
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least a 1— particle cluster (a monomer.) The addition model is not expected to 
support self-similar behaviour in the sense outlined above because the special 
role played by the monomers implies the dynamics gets frozen when one runs 
out of monomers, and the limit state will not be related to any universal pro- 
file (however see [7].) Thus a prerequisite for a meaningful study of similarity 
behaviour in the addition model is the consideration of some mechanism that 
constantly provides new monomers to the system. The most widely studied of 
these mechanisms is fragmentation [2] in which case difficult problems concern- 
ing the large time evolution of large clusters, not dissimilar from the dynamic 
scaling behaviour refered to above, have only recently started to be tackled rig- 
orously |2U1 127) . Another mechanism, is to externally provide for the increase 
of monomers by adding to the right-hand side of the monomer equation in J3J 
a source term Jo(t). A case like this was actually what the original system pro- 
posed by Becker and Doring amounted to ["3EE1, since they considered the case 
where the concentration of monomers stay constant in time. Physically this can 
only be implemented by coupling the system to a monomer bath of infinite size. 
Another way to externally supply monomers, which is much more reasonable 
in a number of pratical applications, such as the modelling of polymerisation 
processes, is to a priori give the source term Jo(£), independently of the existing 
cluster concentrations. The simplest such situation is when Jo(t) = constant, 
and this has indeed been used in applications, for instance in mean-field models 
for epitaxial thin layer growth [3*1 151 1131 115| . 

The study of these Becker-Doring, or Smoluchowski, systems with input 
of monomers, with or without fragmentation, has greatly progressed in the 
mathematical modelling literature (see, for example, [2*3 E3 EU1 02] •) A very 
recent work by one of us [HI] provides a fairly extensive study in the case of 
constant coagulation and fragmentation coefficients and monomer input rates 
given by Jo(t) = at^ . A variety of possible similarity profiles was observed, 
depending on the balance between coagulation and fragmentation and also on 
the rate of monomer input, u>. Our present paper has a much more limited goal: 
we intend to rigorously analyse the constant coagulation (aj — 1) addition model 
with constant monomer input (uj = 0), namely 

' oo 

t Cj = CiCj-i - cicj, j > 2 

where a > is independent of t. 

This is a very special case of the systems studied by [3""], but it is, not only 
sufficiently simple to allow a rigorous mathematical analysis of the similarity 
behaviour to be performed, even including an higher order analysis, but still is 
of some interest for applications [5J. Anyway, it should just be considered as a 
first step towards a rigorous mathematical analysis of the behaviour uncovered 
in the mathematical modelling literature cited above. 

A point that is a clear consequence of our analysis is the way the dynamical 
behaviour of solutions to this infinite dimensional system is really determined by 
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a very small dimensional quantity, namely by the dynamics in a unidimensional 
centre manifold. This reduction of the determining modes of the system is 
the mathematical counterpart to the loss of memory of the initial data implied 
by the physical assumption of convergence to self-similarity. It is reasonable to 
believe that this will also be the case in more complicated systems, but actually, 
at present, this geometric approach to self similarity in coagulation problems is 
not yet entirely clear. 

We now describe the contents of the paper, and briefly discuss their import. 

Our main goal is to understand the large time behaviour of solutions to 
J2J , in particular the precise rates of convergence and possible existence of self- 
similarity. For a solution to exist we require that the sum appearing in @ be 
convergent. In particular, it need not have finite mass M\(t) :— Y^jLi J c j(^)- ^ 
we introduce the total number of clusters as a new macroscopic variable Co(t) 
defined by 



and formally differentiate termwise, we conclude that Co satisfies the evolution 
equation 6q = a — cqC\. Thus, system (0 can, at least formally be written, in 
closed form, as 



In the next section we show that our formal calculations are justified and 
that the solutions of system (0J are in fact solutions to the original system © . 
This is done by the use of a generating function. 

The study of the long time behaviour of solutions requires the knowledge 
of the behaviour of c\ and this can be obtained through a detailed analysis of 
the two-dimensional system for the monomer concentration and total number of 
clusters. This is done in Section [3| and involves the use of invariant regions and 
of a technique related to the Poincare compactification method followed by the 
application of centre manifold theory. The analysis of the full system (0} is made 
possible by the fact that an appropriate change of time t \— > r changes the Cj 
equation to the linear differential equation Cj ' = Cj-i — Cj from where we easily 
get a representation formula for Cj in terms of the non monomeric initial data 
and of ci (see expression JSJ in Section|21) It is this fortunate fact, together with 
the information on c± proved in Section |3 that allow us to study the long time 
behaviour of the other components Cj (t) establishing that all decay as t -1 / 3 when 
t — ► oo , for all j > 1, which will be proved in Section^ The same approach, but 
involving a higher degree of technical difficulties, is used to prove the convergence 
of solutions to a self-similar profile, namely Cj(r) ~ T~ 1 / 2 $i(j/T), and the 
precise determination of the function $i. This is valid for all solutions with 
initial data bounded above by pj~^ with fj, > 1/2, which clearly includes all 
the physically interesting finite mass initial data cases. The precise statement 
and proof takes the whole of Section [S] The main technique involved is a 
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controlled asymptotic evaluation of the sum and the integral appearing in the 
representation formula for Cj. It turns out that the similarity profile < &i(?7) has 
a discontinuity at rj = 1 (see Figure 1.) 




Figure 1: Graph of the similarity profile <I>i (77) . 

To help understand why this is so, it is perhaps interesting to call the reader's 
attention to an heuristic, non-rigorous, but nevertheless enlightening, geomet- 
ric picture of what is happening when we look for similarity behaviour in our 
system. Considering the Cj system for very large j and assuming our level of 
description is such that we can consider j as a continuous positive real variable, 
the equation for c becomes the conservation law 

d T c + djC = (j, t) e M+ x R+. (5) 

Propagation of initial and boundary data for JHJ) occurs along the characteristic 
direction j — r. This implies that initial data for (JjJJ, defined in M + x {0}, is 
not at all felt in the region of M. + x K + with rj := j/r < 1, and reciprocally, 
boundary data for (JSJ, given in {0} x R+, is not felt where rj > 1. Since, in the 
present case, the search for self-similar behaviour entails us looking for limits of 
the solution along lines of constant slope j/r, the geometric picture provided by 
the conservation law leads us to anticipate that the limits of our representation 
formula for Cj, taken with rj < 1 will demand less stringent conditions on the 
initial data than limits with rj > 1, and reciprocally for the requirements on 
the "boundary data" c,\ (r) . And of course the border case rj = 1 is expected to 
exhibit some problems. This turns out to be exactly what happens. 

Naturally, the discontinuity in <f>i is an indication that along the character- 
istic direction rj = 1 the solutions do not scale as t~ x / 2 (if they scale at all.) In 
fact, we end the paper by proving, in Sectional that there is another similarity 
variable, £ := ^7=-, that essentially corresponds to a kind of inner expansion of 
the characteristic direction 77 = 1, and a similarity profile, $2(^)1 such that, at 
least for monomeric initial data, solutions behave like Cj{r) ~ t _1 / 4 $2( :2 -^"), 
where the function $2 (to be defined in Section^ can be expressed jHJ m terms 
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of Kummer's hypergeometric functions P pp 504-515], and have the graph pre- 
sented in Figure 2. It is worth noting the similitude of the graphs of the scaling 
functions in Figures 1 and 2 to those in [SJ Fig. 2] and Fig. 2, 8]. 




-10 -7.5 -5 -2.5 2.5 5 

Figure 2: Graph of the similarity profile 3>2(£) 

The restriction to monomeric initial data is probably not important, but we 
could so far not overcome the difficulties involved in controlling the initial data 
term for small cluster sizes. We discuss these difficulties at the end of Section|SJ 

Putting together the results in the last two Sections of this paper, we can 
conclude that solutions do converge to a similarity profile, that is a unimodal 
distribution with a spike located at j ~ t 2 / 3 decaying to zero like i -1 / 3 for js 
behind the spike, and the spike itself decaying like i -1 / 6 . This means that the 
family of solutions {cj(t)} tends to become more spiky as we look for larger and 
larger times t and cluster sizes j. 



2 Equivalence of the related system 

Let (cj(i))°^ be a nonnegative solution of Q. To show that this is also a 
nonnegative solution of the original addition system requires that we prove 
the convergence of Y] , c n (t) to Co(t). To this end, it is convenient to introduce 
a new time scale 

r{t) := f Cl {s)ds (6) 
Jo 

along with scaled variables 

S,-(r):=c,(t(r)), (7) 



G 



where t(r) is the inverse function of r(t). Since c\{t) > 0, these are well defined 
and t is an increasing function of t. Now consider the Cj -equations in flljl: 

6j = CiCj-i - CiCj, j > 2. 

By the change of variables (t, Cj(t)) i— > (r, Cj(r)) these equations can be written 
as a linear system 

V = Cj-i -cj, j>2, 

where (•)' = Now this system of ordinary differential equations is a lower tri- 
angular system that can be solved recursively starting with the j — 2 equation. 
Doing this we obtain, by the variation of constants formula, 

c,(r) = e- £ jf^CkiO) + JJ 1 ^ jf Mr - s^e-da. (8) 
We now prove the following. 

Theorem 1 If the series of initial data is convergent, ie if Yl'jLi c j < °°> 
then a solution of system ^ will also be a solution of the original system J3J). 

Proof 

Introduce the following generating function: 



71 = 2 



F(T,z):=J2cn(r)z n . (9) 

n=l 

Using (JSJ) we can rewrite this as follows: 

oo 

F(t, z) = c 1 {T)z + ^2c n (T)z n = c 1 (t)z + G(t,z)+H(t,z), 

where 

G(r, z) := e~ T ^ l__f_ Cfc (o), 

7i=2fc=2 v ' 
00 „n /-T 

n=2 [ >■ Jo 

Examining these expressions, we find that the series can actually be summed. 
For G we have 

00 n ~ 2 n-k-2 n 00 00 m m+k+2 

G (r, Z) = e - r EE ^^Mo) = e ~ r EE — ^r-^(°) 

71=2 fc=0 v y fc=0 m=0 



°° " T n-k z ,. 



T E ( E ^- ) wo) - e - T(1 - 2) E ***(°)- 

fc=0 \rn=0 ' / fc=2 
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By hypothesis the series in the above expression converges for \z\ < 1, so we 
have 

G(t,z) = e- T ( 1 " z )(F(0,z)-c 1 (0)z) for \z\ < 1. (10) 
As for H we have 

= z 2 f T c 1 {T-s)e sz e- s ds = z 2 f Z^e-^-^-^ds. (11) 
Jo Jo 

The expression for F now becomes: 

F(r, z) = c x {t)z + e"^ 1 "^ (F(0, z) - c x (0)z) + z 2 f c 1 (s)e-^~ s){1 - z '> ds, 

Jo 

(12) 

which, at z = 1, yields 

F(r,l) = c 1 (r)+F(0,l)-c 1 (0)+ / ci(«)tk. (13) 

Jo 

One can easily verify that F(t, 1) given in (|13f) satisfies the first equation of JIJ 
which proves that F(t, 1) = c (r). I 

Having established that solutions of J2J are in fact solutions of ©, all sub- 
sequent analysis will be carried out on system Q. 

3 The bidimensional ODE system governing the 
monomer dynamics 

From the reduced system (0J we observe that the equations governing both the 
monomer dynamics and the total number of clusters are actually the bidimen- 
sional (co, ci)-system, the dynamics of which can be studied quite independently 
of the remaining components. To simplify notation we will use x = c\ and y = cq 
in the study of this bidimensional system. 

Let a > be a constant, and consider the system 

V = a-xy 
x = a — x — xy 

We are interested in nonnegative solutions to (|14f> and so, from hereon, everytimc 
we speak of solutions we actually mean nonnegative solutions. Our first result 
concerns the gross features of the long time behaviour of solutions. Another 
result, to be presented in Proposition |2 will establish some finer details of the 
long time behaviour. 
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Proposition 1 For any solution (x, y) of {1$ the following holds true as t — > 
+00: x(t) — * 0, y(t) — ► +00, and x(t)y(t) — ► a. 

Proof: Let fl be the connected subset of M. + x M+ whose boundary is {y — 
0} U {x = 0} U {xy = a}. Consider first our initial data in ft. Elementary phase 
plane analysis shows this set is positively invariant for the flow (see Figure 3.) 




Figure 3: Region O with a sketch of the flow in 9fi. 

We now observe that 1|14[) does not have equilibria, and consider the subset 
0,1 of fl defined by 

ill := \(x, y) G O : max |o, x\ < y < — \ . 

v v 20 J 2C X 

Again, elementary phase plane analysis shows that Cl\ is positively invariant 
and, for any initial condition in £1 the corresponding orbit will eventually enter 
f2i (see Figure 4.) 

From this we immediately conclude that, as t — > +00, we have x(t) — > and 
y(t) — * +00. Furthermore, for all initial data in O, there exists a T (depending 
on the initial condition) such that, for all t > T, the orbit is in fix, and so 

a a 2 

t > T =>• x < y < — <^>a — x < xy < a. 

x x 

Since we know that x(t) — -> as t — > +00, applying limits in the above inequality 
gives 

lim x(t)y(t) = a. 

t — >+oo 

Consider now initial data (xo,yo) E Cl% = R + x ]R + \ fi. Fix K\ > Xq, 
K 2 > y and let ^2(^1,^2) = ^2 n ([0,ifi] x [0, -Kjj])- % the analysis of the 
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xy = CX-x 2 



Figure 4: Region ill with a sketch of the flow in dft, dfli and in il. 



flow in 9^2(^1,^2) we conclude that the orbit will eventually enter f2i (see 
Figure 5) and so the previous analysis apply. 

This concludes the proof. I 

In the next Proposition we prove results about the rate of convergence of 
x(t) and y(t) as t — ► +00, using an approach akin to Poincare compactification. 

Proposition 2 For any solution (x, y) of \1$ we have: 

/3 \ 1/3 

(i) lim -t x(t) = 1 
t->+oo \a J 

(ii) lim (3a 2 ty 1/3 y(t) = 1 

'+00 



3 ^m 



(iii) lim -t (a - x(t)y(i)) = 1 

Proof: For the present proof we found convenient to work with the variable 
v := a — xy instead of y, and with a new time scale. The idea is essentially that 
used in Poincare compactification: to bring to the origin the critical point at 
infinity, and to turn it into a point of the phase plane by appropriately changing 
the time scale (cf. eg O chap. 5].) By changing variables (x,y) 1— > (x,v) 
system (|14t becomes 

' x = v-x 2 

v = ax-2xv-a^ + ^. 1 ' 
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Figure 
O. 



e 5: Flow in the set &,2(Ki, K2) surrounding an initial point (xq, yo) outside 



Suppose x(0) 0. Otherwise, since we know by the phase plane analysis in the 
proof of Proposition^that x[t) > for all t > 0, just redefine time so that x(0) 
becomes positive. Now change the time scale 



x(s) 



ds. 



and define (x((),v(()) := (x(t(()) 7 v(t(())) , where t(£) is the inverse function 
of C(t)- By Proposition ^ we have x(t) — > as t — > +00, and so also £ — ► +00 
as t — > +00. With the new time scale system l|15ll becomes 



X 


1 


' 


" 




X 


+ 


V 







—a 




V 





~2 1 ~2 

ax + v 



2vx z 



(16) 



where (•)' = Jr. The region of interest, corresponding to (x,y) G K + x R + , is 

(x,v) € R + x (— 00, a), but actually (|16fl is valid in all of the phase plane K 2 . 
In this way, the critical point at infinity of system l|14|) is mapped to the critical 
point at the origin for i|ltj|) • From the results in Proposition ^ we know that 
all orbits of (|l(jfl obtained from orbits of (|14fl by the above map will eventually 
enter M + x (0, a) for sufficiently large times £ and converge to (0, 0) as C — * +00. 
Using standard results chap. 2] it is straightforward to conclude the existence 
of a centre manifold for (|16|l given by 



v = 4>{x) where (j)(x) = x 2 x 4 + 0(x 6 ), (17) 

a 

that locally attracts all orbits. The dynamics on the centre manifold is given by 

0(x 7 c ) as C -> +00. 



a 
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Write this differential equation as 
d f a 



1 + 0(x 2 c ) as C -> +00. (18) 



Fix an initial condition (Coj^o) and fix e > arbitrarily. Since we know that 
^c(C) as £ — > +00 we conclude that there exists a T > Co such that, for all 
C > T, the following inequalities hold 

_d_ / q \ 



Integrating these differential inequalities between T and £ we get 

+ (l - e )(C - T) < JL < Jr + (l + £ )(C - T), (19) 

where 2?t = ^(T 1 )- Dividing (|19fl by £ and taking lim and lim we obtain 

a a 

1 — £ < hm — < lim —^7 < 1 + e 
~ c~oo 4C^ " C-+00 ACxt ~ 

and since e > is arbitrary, we conclude that 

/4 \ 1/4 

hm -C i c (0 = l (20) 

for solutions corresponding to orbits of ill6|) on its centre manifold. From stan- 
dard centre manifold theory .8., chap. 2], long time behaviour of (x((),v(()) 
is determined by the behaviour on the centre manifold modulo exponentially 
decaying terms O (e _A ^) , where A € (0, a), in particular we can write 

Multiplying this equality by (^C) 1 ^ 4 an d using (|2U|) we conclude that 

/4 \ 1/4 

lim -C £(C) = 1. (21) 
C^+oo \ a J 

We can do the same for the evolution of the v variable: for orbits on the centre 

1/2 

manifold we have v c — x 2 c + 0(x*), multiplying this by (— C) 1 an d using l|2U|l 

we conclude that 

lim -C « C C) = 1. 

To obtain the dynamics outside the centre manifold we proceed as with the x 
variable and get 

4 \V2 



hm -C «(C) = 1- (22) 
C— »+oo V a 
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In order to obtain the corresponding rate estimates in the original time variable 
t we need to relate the asymptotics of both time scales. By the definition of the 
new time scale 

C(t) - C(*o) = / ~^ds 

Jto X{S) 

and so ^ = j 1 ^-. Hence, by the inverse function theorem, ^ = x((), and thus 

*(C) - *(Co) = f x{s)ds (23) 
J Co 

where Co = C(^o)- From (|21() we know that 

Ve > 0, 3T = T{e) : VC > T, f ~cj x(£) G [1 - e, 1 + e] (24) 

Let us look first at the upper bound: for all £ > T, 

/4 V 1/4 
2(0<(rC + ^ 



substituting this into l|2"3*)) we get, for (o > T, 

i-C / 4 X-l/4 

*(C)-*(Co) < y (-*) (l + e)da 

Multiplying this inequality by | (j^^C -3 ^ 4 an d passing to the limit as £ — > 
+oo results in 

Em |(i) 1/4 r ^(C)<l+ £ . 

Using the lower bound of (|24|) and the same argument we obtain the reverse 
bound 

Hm |(iy /4 c - 3/4 *(c)>i- £ . 

and by the arbitrariness of e > we get 
3 /4\ 1/4 

C" 3/4 i(C) = 1 + o{l) as C -» +oo, (25) 



4 \a 

which, together with if^Tl) and allow us to conclude that, as t — > +oo (ie, 
as C — > +oo) 

3 \ 1/3 / 4 X 1/4 



t x(t)=[-C) i(C)(l + o(l)) = l + o(l), 
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g x 2/3 / 4 x 1/2 



and 

' ^) w( * )= t?) «(C)(i + c>(i>> = :i+^(ii. 

thus proving (i) and (iii) respectively. 

Finally, to prove (ii), observe that for the original variable y(t) we have 
y = 9 ^r E and thus, as t — > +00, 



(if) 1 ' 3 -!*) (f<) 1/ -(*)a(f') 2/3 
which concludes the proof. 



4 Long time behaviour of the system 

We now turn our attention to the long-time behaviour of Cj(t) with j > 2. It is 
convenient to consider the time scale introduced in Section [21 The first result 
we need is the following, the proof of which is entirely analogous to that of the 
corresponding results presented in the second half of the proof of Proposition^ 
and so we refrain from repeating it here, 

Proposition 3 With (cj), r, and (cj) as given in and Q), the following 
holds true: 

(i) lim \ (-X r a /3 r (i) = 1; 

fir ~ 

(ii) lim \ — ci(t) = 1; 
t^+oo V a 

With this knowledge, we can now use (JSJ) to get information on the long 
time behaviour of Cj(r) for j > 2, specifically, we prove now that (ii) in 

Proposition holds for all Cj(r). To see this, multiply © by y^^- The first 
term in the right hand side of the resulting expression is the contribution due 
to the non monomeric initial data, and since j is fixed, we have, as r — ► +oo, 

fc=2 U 

for every A < 1. In order to study the remaining integral term, first change the 
integration variable s h- > y — s/t in order to get the integral in a fixed bounded 



region, and define the function ip(-) := y ^ (•) c i(0- We thus get 
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Let < e < 1 be fixed, and write the integral as J Q £ + f l e . The last integral is 
fairly easy to handle: Since ip(s) is a continuous function and, by Proposition^ 
is 1 + e>(l) as s ^ +oo, we conclude that it is bounded in [0,oo) and so there 
exists a positive constant such that < ip(s) < for all s. From this it 
follows that 

i-e vi-j/ JisVi-y 

< M^e^ 1 -^ f ~^=dy 

Ji-e V 1 - y 

and so it is also exponentially small when r — > +oo. For the integral over (0, 1— e) 
we use y < 1 — e => r(l — y) > re — * +oo as r — > +oo, and Proposition^ (ii) to 
conclude that ip = l + o(l) in the region of integration, provided t is sufficiently 
large. Therefore, V<5 > 0, 3T(S) : Vr > T(5), ip(r(l - y)) € [1 - 5, 1 + 5], and, 
as t — > +oo, 



1 " e V>(t( 1 - y))y ^ 



(1 - S)I 0J (r) < / nH ™ e-^dy < (1 + J)/ 0ji (r) (26) 



with ^ 

I 0j -(t):= / "^Le-^dt/. 

,/o vi - y 

This integral is easily estimated using Watson's Lemma [J, pp 427-8]: from the 
Taylor expansion 

y>- 2 ,- 2 y^ (2fc-l)!! k 

which converges for \y\ < 1, we can apply Watson's Lemma directly to get 

Ioj(t) = r(j ,I 1 1) asr^+oo. 
Plugging this into (J25J) immediately results in 

/ 1-£ ^(1-10)^-%^ = x + (t - 1} asr ^ +oc . 

Jo V 1 -y 



(i-2) 

This, together with the other two exponentially small contributions proved 
above, implies that Cj-(t) — > 1, as r — > +oo, for all j > 2, thus com- 

plementing part (ii) of Proposition [3] Using part (i) we can easily translate 
this behaviour to the original time scale, and conclude the following: 

Theorem 2 Let (cj) be any non-negative solution of |3]) with initial data sat- 
isfying cq(0) = ^2j—i Cj(0) < oo. Then, as t — > +oo, we have 
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/3 \ 1/3 

(i) -t Cj (t) — >1 for all j > I; 



a 



(ii) (3^r 1/3 $>(t)^l; 




5 Self-similar behaviour of the coagulation sys- 
tem outside the characteristic direction 

We have now reached the main objective of the paper: the study of convergence 
to similarity profiles in this addition model. 

Let $1 : R+ \ {1} — > R be the function given by 

1 

if V < 1, 



if 77 > 1 

The main result in this part is the following: 

Theorem 3 Let (cj) fee any non-negative solution of |3J) with initial data sat- 
l 

2 

and 0) respectively. Then, 



isfying Bp > 0, /i > ^ : Vj, Cj(0) < p/j M . Let r(t) and Cj(r) fee as given in 0) 



/ 2 _ 
lim a —t Cj(r) — 

rj—j/r fixed 

Proof 

Expression (jSJ) describes the time evolution of c,-(r) by the sum of a contri- 
bution dependent only on the non-monomeric part of the initial data with a 
term determined by the behaviour of C\{t) in the appropriate time scale. For 
monomeric initial data, only this last term is relevant, and we start our analysis 
by it. 

5.1 Monomeric initial data 

Assuming monomeric initial conditions we have, for j > 2, 



2 ^ 
TMT) =U-2) 



I c 1 (T-s)s j - 2 e~ s ds. (27) 
Jo 



1G 



Consider the function ipi defined in [2, oo) x [0, oo) by 



^{x,r):=-f -/ c 1 {r- S )s x - 2 e- s ds. (28) 

r (z- 1) Jo 



Clearly ipi — J c if 2 < x = j e N, and the use of tpi is a good deal 
more convenient in order to obtain the required asymptotic results. So we now 
concentrate on Let r\ :— x/t. By changing variable s i— ► y = s/t, using the 
recursive relation T(x — 1) = T(x)/(x — 1), and Stirling's asymptotic formula 
T(x) = e~ x x x ~z v2tt (l + 0(a;~ 1 )) as a; — * oo, we can write, as r — * +oo, 

yi(*?r,r) = j= (l + Oir- 1 )) j y -^— 2 dy. 



V 

J29) 

In order to make use of our knowledge about the long time behaviour of c\ (r) 
given in Proposition 3 it is convenient to rearrange the integrand of l|29[l by 



multiplying and dividing it by y ~(1 — y), resulting in the following, as r — > 

+oo, 

Vl (r,T,T) = ^=r,z-^Tz{l + 0{T- 1 )) / ^(r(l-y)) 2 dy, (30) 

v ^tt jo y v 1 — y 



where -0(0 := y „ (') c i(') i s the function that was already defined and used in 
Section @] The proof now reduces to the asymptotic evaluation of the function 

I(r h t) := r-*ir nT e Tri / 0(r(l - y))-— = dy (31) 

jo y v 1 — v 

as r — > +oo. 

We first consider the case 77 > 1. 

Observe that y-SgT^logy-y) _ ^r^logy-ry = . ^(y^ where ^ . R + _^ R ^ 



2 

7,-1 



defined by the last equality, satisfies, for all y € (0, 1] and r > 

9'x{y) = - r > (rry - 2) - r > 0. 

2/ 

Therefore, 31(3/) < <7i(l) = — r, and so 

e T(-q\ogy-y) pi i 

^(r(l - y)) dy < Afye" T / ^^dy = 2M^-. 

y 2 V 1 ~y Jo -y 

Plugging into J^U w e obtain, as r -> +00, /for) < CV^e^^ 10 ^^^ for 
some positive constant C. Since 77 log 17 — r\ + 1 > for all rj > 1, we conclude 
that I(rj, t) — ► as t — ► +00, and this proves the result when 77 > 1. 
Consider now the case r\ £ (0, 1). 
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Observe first that the exponential term inside the integral in l|31l) has a 
unique maximum at y = rj. Therefore, in order to estimate l|31[) when r — ► +00 
we are going to study its behaviour around the maximum separately from the 
other cases. To this end, let e < min{?7e , 1 — 77} and write 

7(ry,r) = r-rf^e^ (J +J +J J ^(r(l - v ))-_==-dy 

=: h(ri,T)+I 2 (r,,T) + h(ri,T). (32) 

The study of the integral h(r), t) is entirely analogous to what was done in the 
case T] > 1 : since < y < e < rje^ 1 < 77 we now have, for all r > jjz^tj^, 

, r/T-2 t{t] - y) - 2 

ffi(y) = r = > 0- 

Thus, gi(y) < ffi(e) < gi(v e l ) — ( 7 7 r — 2) log 77 — (777- — 2) — rr/e^ 1 , and we can 
write 

Jo Vi-v 

< CT^e- Tri/e — >0 asr^+oo, (33) 

for some positive constant C. 

For the integral 73(77,7-), define the function g 3 : (0,1) — > R by ga(y) := 
(77 log 77 — rj) — (77 logy — y). Since 53 has a unique minimum at y = 77, where its 
value is 33(77) = 0, we conclude that y > 1 — e > 77 => #3(7/) > y 3 (l — e) > and 
thus, there exists a positive constant C such that 



e --rga(,y) 

h(v,T~) = r* / t/>(t(1 -y)) 2 ^ dy 



yV 1 - y 



< CT^ e - r93(1 - £) — >0 asr^+oo, (34) 

From l)33|) and (|34|l we conclude that 

1(77, t) = 72(77,1-) +o(l) asr^+oo. (35) 

In order to estimate 12(77, t) we proceed as follows: we now have r(l — y) > 
re — > +00 as t — > +00, and thus iP{t(1 — y)) — 1 + e>(l) when evaluating 7 2 for 
large values of r. As in a similar situation in Section^] we have V<5 > 0, 3T(5) : 
Vt > T((5), t/'(t(1 — y)) e [1 — 5, 1 + S], and we can estimate 

(1 - (TIt^-^ J(?7,t) < h{r],r) < (1 + S)T^T 1 - rlT J(T 1 ,T) asr^+oo, (36) 
with 

J(77,r) := / 9 2 e- T ^y) d|/> 

and ^ : (0, 1) — > R is defined by 0(y) := y — 77 logy — 77. Since this function 
is smooth and has a unique minimum, attained at y = 77 G (e, 1 — e) with 
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value <j>{q) — —77 log 77 and 4>"{rj) = J/ --1 , Laplace's method for the asymptotic 
evaluation of integrals |5J pg 431] is applicable to J(rj,r), and we obtain, as 

T — ► +OO, 

1 / 2-7T / P TV log 1 \ 

J{ "- r}=c '" "wmivj7, +0 {— )■ < 37 > 

Hence, from (jSUJ), JSSJl, and J57JI, we can write, as r — > +00, 

Mnr,r) = -L^rV^ 2 I J ^- (l + O (r" 1 )) (38) 

+ ^ r '"" rT5 ° (-^75- J ( x + ( r )) + °w ( 39 ) 

After a few trivial calculations, we imediately recognize that, as r — > +00, 
(H) is equal to ^^=(1 + o(l)), and ® is r^Ofl) (l + (t" 1 )) + o(l) = 
O (t- 1 ) + o(l) = o(l). Therefore, for 77 G (0, 1), 

ipi [rjT, t) — > as t ► +00 

and this concludes the proof in the case of monomeric initial data. 



5.2 Non monomeric initial data 

If the initial data has non zero components Cfe(0) with k > 2, the proof of the 
stated similarity behaviour ofcj(r) requires, according to (JHJ), that we now prove 
that 

[2 _jL, T J—k 

. lim J -r e- T 7— MT c fe(°) = °- 

J,T^+OC V ol [j — k)\ 

rj—j/r fixed k—2 

Define v := r/^ 1 , write r = jv, and use the assumption on the initial condition, 
namely Cfe(0) < p/k^. We thus have 



n L — ' (7 — fc)! \ a (7 — k)\k^ 

k=2 w 7 fc=2 w ; 

a 

where <y9 2 is defined by the equality. Our goal is to prove that <^ 2 (i / , j) — > as 
j — > 00, for all positive 1/ ^ 1. 

Like it was done before, we shall study the cases v > 1 and t/ < 1 separately. 
By the heuristic geometric reasons explained in the Introduction, the case v > 1 
(77 < 1) is now easier to handle. In fact, for this case, even using the notoriously 
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bad upper bound to the sum provided by taking the product of the number 
of terms by their maximum, we only need to impose \i > 0. Using this same 
approach to estimate if 2 in the case v < 1 (rj > 1) we would need to impose 
fi > 1, which would not even guarantee that all finite mass initial data are 
included. So, for this second part of the proof, a finer approach is required, which 
consists in estimating the small and the large £ contributions to if2 separately. 
Consider first the case v > 1. 

Change the summation variable k t— > £ :— j — k. It is sufficient, for this range 
of v, to bound i/?2 as follows 

9><y, J) = ^ E < e-» E (40) 

( ■ \* 

Considering the sequence ui := ^ , and studying the sign of 

we conclude the maximum of ui is attained at I — \jv\ > jv— 1 > j — 1 > j — 2. 
As j — > 00, we can thus estimate the right-hand side of 140J1 : 

£=0 



where we use Stirling's approximation in the second equality. Since v > 1, we 
have 1/ — 1 — log > 0, and so —j(y — 1 — log i/) + log j — > —00 as j — > 00 from 
which we conclude l|41l) goes to zero as j ' ^ 00, and we obtain the result we seek 
in the case v > I. 

Consider now the case v £ (0, 1). 

Let (3 € (ve x ~ v , minjVe, 1}) be fixed, and write 

\0<.6</3? u 7 0j<t<j-2 KJ > J 



< — V^e 



'(i-i) 



ii- 2 



(i-2)! 



t/- 3 /a 
2' 1 v / 27r 



-3/2 -j* 



(j - I)' 



j 3/a e -i-y_i) 



C? - 1)1 

2 W 



2^V2tt 



i-i 



e ^j-^j ^ 

(41) 
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For the first sum we have 



jv e 



o<e<Pj w yjJ 



E 



J v r > o<e<pj 



< \fjv e 3 V 1 



/77 i _ l 
-j 2 p — ► as j ^ oo, if /x > -. 



By Stirling's expansion we can write, for all sufficiently large t, l\ > e~ e £ e+ ^. 
Therefore, since in S2 we have (j — > 2 P > 1, we can estimate that term as 
follows, 



< 



/3j<e<j-2 e 4 



~ -> ST o-k ( 3 ve 

0j<e<j-2 

V" S (x 

/3j<e<j-2 v 

, ,i-2-L«J 
1- tt /„„\ LA?J+i 



v -jv \ P I ive 



/? 1-f 



^ /? (ve 



A-v \ 3 



ve — (3 ve \ /3 
This completes the proof of the theorem 



as j — > oo. 



6 Self-similar behaviour of the coagulation sys- 
tem along the characteristic direction 



Consider the function $2 : K — > K defined by 



$ 2 (£) : = e -i? 2 / e-^-^dw. 
Jo 



The main result of this part of the paper, complementing Theorem [3J is the 
following 
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6.1 Monomeric initial data 

Theorem 4 Let (cj) be any non-negative solution of @) with monomeric initial 
data. Let r(t) andcj(r) be as given in ^Bjj and respectively. Then, 

\ 1/4 

. lim — r c i (r) = $ 2 0- 

£=i=5 fl xo d 



After a few manipulations with the integral defining $2 we can write 

= J^^ ?2 (J-* (iC 2 ) +sgn(-0/x (i£ 2 )) , 

where /„(•) is the modified Bessel function 1 , pp 374-7], and the value at £ = 
is defined to be the limit of the right-hand side as £ — * 0. Actually, $2 can be 
seen as an element of a larger family of functions: it is the function obtained by 
making v = | in 

■■= 

_ T(l-^m\e-^ 2 ( , _ 1 3. 1 2) sgn(-Q + l M ft _ 1 3. 

These functions were formally deduced in |31j as similarity profiles for addition 
models with time dependent monomer inputs Jq{£) = at u , where v — — 
The functions £7 and M are Kummer's hypergeometric functions 1 , pp 504-515]. 

Proof 

For monomeric initial data JSJ reduces to 

7T 2 \ 1/4 ( -"2"T J r r 

Consider the function ip^ defined in [2, 00) x [0, 00) by 

1/4 



V3(z,t) := YT^Ty / - s)^ 2 e- s ds. (42) 

In order to consider the similarity limit we introduce the variable £ := and 
write l|42(l in the similarity variable £ as follows: 



I(t + £V t - 1) Jo 



(43) 
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( 2 \ 1/4 ~ 

Since </?3 (j, r) = ( jjrr) Cj (t) if 2 < a; = j 6 N, to prove the theorem it suffices 
to show that 

lim p 3 (t + £Vt,7-)=* 2 (0- ( 44 ) 

r — >+oo 

We start by considering, in the integral in the right-hand side of l|43|) . the change 
of variables 

\ T 



and by writing if>(-) := J j^(-) ci(-). This results in 



,1/4 



1 (r + - !J Jo 

From Stirling's expansion for the Gamma function we have, as t — > +oo, 
e-^/^T + fv^ T+ ^3 Jo V V^/ 

(45) 

We can now estimate the behaviour as r — > +oo of the multiplicative factor in 
the right-hand side of l|45|l as follows: 



f f 



1 TTT- 1 — T(l + O(l)) 



1 + 



e-« 2 (l + o(l)) 



= e5«V« 2 (f + e>(l)) = e"5? 2 (i + o(l)) 

where the last but one equality is obtained by changing variable t^i:= 
and applying L'Hopital's rule twice. Thus, we can write l|45(l . in the r — > +oo 
limit, as 

T 1/4 

^ 3 (r + ^,r)-e-^ 2 (l + o(f)) / ^(^ 2 )e (T+ ^ 2) +V ^^- 



(46) 

In order to estimate the integral, and considering that the only relevant informa- 
tion available about tp(s) is that it is a non-negative, continuous, and bounded 
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function on [0, oo) and tp(s) = 1 + o(l) as s ^ +00, we are forced to treat 
the region with w close to zero separately. The idea is to write the integral as 

T l/4 

j Q + and to prove the first integral can be made arbitrarily small. More 

precisely, we now prove the following: 

V 5>0 ,3 eo=eo(0 :V £<£0 , [ E ^w 2 )e {T+ ^- 2)l < 1 -£) + ^d w< S. (47) 

ro=r (C,eo) T>Ta J 

In fact, let S > be arbitrary, let eq = £o(£) := max|e : e^'^ 2 — < o| , 

with M$ = sup ip(s), and let r := maxMe^, |(VC^+8 - £) 2 }- Then, for 

sS[0,oo) ' 

all t > t we have r + £y/r — 2 > 0, and, using log(l — x) < —x, we conclude 
that, for e < So, 



V(^ 2 )e (T+Sv/F - 2)los ( 1 -^) + ^ 3 ^ < 
Jo 



5 



_l/4 

We are now left with the contribution f to the integral in (|46|l . Since 
to > e y/rw 2 > \fre 2 — ► +00 as r — > +00, it is easy to see that, for all 
sufficiently large r we have the integral of interest asymptotically equal to 

{l + o{\)) f^e^^^-^^^dw. (48) 
To estimate (|48|1 observe that, as r — > +00, 

e (r+Cv^-2)log(l-^i)+^ 2 = 
^ ? / ?/I 2 \ T I / 

1 - 




(1 + 0(1)) (49) 

where the last equality is obtained by changing the variable r 1— > a; := and 
applying L'Hdpital's rule. From this we conclude that there exists a continuous 
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function g(w; r) denned for r > e 4 , w 6 [e, t 1 / 4 ), satisfying 1 + g(w; t) > and 
<?(«;; r) — ► as r — ► +oo for each fixed w, such that 

e (r+^2) lo g (l--|) + V W = e -^-|^ (1 + g{w . r)) _ (5Q) 

Considering r sufficiently large, for instance larger than tq previously defined, 
we can use the bound on the logarithm to obtain, from J5UJ), 

1 + g(w,T) = 

= exp (V + £Vt - 2) log ^1 - + yftw 2 + £w 2 + ^w 4 

// , r- „n / it; 4 u> 6 \ _ , „ , f & 

< exp ^(t + ^ 2) (--^ - — - ^ J + VtW + £™ 2 + 

f n w 2 £ w 4 w A 1 , _ w 6 

= exp 2^-^ + o(t + £Vt-2) 



ft 2 \ft t 3 V sv V/ 2 



_. e Hw;r) 



where the last equality defines h. Change variables (w,t) i— ► (u, 6*) where w := ^= 

and 9 = -^/r. Denote by h the function /i in the new variables. We immediately 
conclude that h is the polynomial 

h(v; 6) = 2v- |u 2 6> + v 2 - \{9 2 + £9 - 2)v 3 

and the region of interest in the new variables is A := {(v, 9) : 9 > e 2 , %■ < v < 
1}. Clearly h is bounded from above in A. Denote by logM an upper bound 
for h. Thus, for all t > To we can write 1 + g(w; r) < M. Since g(w; r) — ► as 
r — > +oo for each fixed w, the quantity 

t w := max{r : ^(w; r)| > 5} 

is well defined for each fixed w and, because g is continuous, 

T-i :=max{ 

} < +°°- 

Two further auxiliary estimates, both immediate, are needed: firstly, knowing 
that e-^-l^ g L 1 (M), we can define 

<>+oo 

l(£,w ) := / e -«- 2 "^ 



and conclude that 

j(e,wo)<<y, (5i) 



25 



and secondly, given 6 > arbitrarily, and defining e\ = £i(£) by the expression 
E\ := max |e : e'^' £ — ~ < o| , it is obvious that, 

Va> 3 ei :V E<ei , / e"^ 'dw < 5. (52) 
Jo 

From the above estimates, and using the notation already introduced, we deduce 
that: 

V5>0)3e=min{e ,ei} ■ ^e<£i ^w-max{£,Wd} • ^wo>wi ^n— ti(to,Wq,e) • ^t>ti j 

TO = To(i,So) 



-1/4 

(1 +g{w 1 T))e- ( ' w2 -^ wi dw 



(1 + g(w; r))e~Z w dw + (1 + g(w; r))e^ w dw 

J Wq 

pWq 

< (1 + 5) e -^ 2 -^ 4 dw + Ml(i,w ) 



r+oo 

< / e-t w2 -? wi dw + (l(Z,0)+M)6 
Jo 

r+ca 

< / e-^-^dw + CS 
Jo 

and also, using the positivity of 1 + g(w;r), 



T l/4 

(l+g(w,T))e-t w2 -? w4 dw > 
> / (l+g(w;T))e-Z w2 -? wi dw 



> (1-8) e~t w dw 



-foo / pe poo fWo\ 

e-e^-Wdw- I / + / +6 )e-^ 2 -i wi dw 

r+oo 

> / e-t w2 -^ w4 dw-(l(£,0) + 2)6 
Jo 

r+oo 

> / e -& 2 -k™ dw- C6 
Jo 

where C := !(£_, 0) + max{2, M}. 

Plugging this result into (|48() and remembering (|47|l we finally conclude (|44|) . 
This concludes the proof. I 



2<i 



6.2 Remarks concerning non-monomeric initial data 

It is our conviction that Theorem 21 also holds for non-monomeric initial data 
satisfying the decay condition of Theorem A few numerical runs gave results 
consistent with this conviction. However, we were not yet able to prove this is 
indeed so. In the present final section we shall describe what was achieved so 
far, and what were the difficulties encountered. 

From expression © , if we consider non-monomeric initial data in Theorem^ 
we need to study the value of the limit 

&) 1 '\-± T f^ cm . (53) 

£=i=t fixed k = 2 

Using the definition of £ we can write j = t + £\/t, which, solved for r gives 
t = jAj where the function Aj : R — > R is defined by 





A i(0 := h 1 + 77 +sgn(-0J7- ■ (54) 



A number of properties of this function can be easily established, in particular, 
f < 1 if £ > 

Aj < = 1 if f = (55) 
[ > 1 if £ <0, 

Aj — ► 1 as j — > oo, for all £, (56) 
and the following power series expansion 



A - = 1 —f + —f 2 + V( i r -i (2"-3)!! (C\ - for , > V 



(57) 



Assuming Cfe(O) < p/k^ , changing the summation variable to £ := j— k, writing r 
in terms of j as indicated, and not taking multiplication constants into account, 
the limit (|53|l reduces to 

lim jVVjA, y JiM. . (58) 

which corresponds to the study of tp2 in the proof of Theorem [21 The fact 
that now Aj is not bounded away from 1, as was the case with v in the above 
mentioned proof, is the origin of the main difficulties. The most promising 
approach to the evaluation of (|58l is to proceed by decomposing the sum into a 
small I and a large I contribution, similar to what was done in the case v € (0, 1). 
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In that occasion the cut off size separating the two sums was at I = (3j for a /3 
fixed in a given way. This cut off scales like r as a function of j. It is natural 
to keep considering a cut off fulfilling the same scaling requirement, which in 
this case means ~ j — |£|VJ- So let us define j* := (j — 2) — (1 + |£|)VJ — 2 and 
write the expression in 1)58(1 as 

\0<i<j* w y j*<i<j-2 w ' / 

The small ^ sum, £3, that corresponds to the contribution of /arge cluster in 
the initial data (remember the change of variable k t—> £), can be estimated in 
the same way as was done for the sum Si in the proof of Theorem [3] 

S*(i) < fV^e-iAi V" OA)* 

0<*<j* 



l/2\^ # 



E 



e 

^ + (i + iei)(i-fj 1/ y °<<<i* 

< QHe^X^.QH^o as^oo, if M >l 
f=o 

The large € sum, 54, corresponding to the contribution of small clusters, is 
the one that could not yet be tackled rigorously and we must at present leave 
it as a 

Conjecture 1 With the definitions above it holds that S^i^j) + S^j) — > as 
j — > oo, u>/ien /i > i, and hence Theorem^also holds for non-monomeric initial 
data with this decay. 
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